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ABSTRACT 



Recent studies suggest that infrared dark clouds (IRDCs) have the potential of harboring the earliest stages of massive star formation 
and indeed evidence for this is found toward distinct regions within them. We present a study with the Plateau de Bure Interferometer 
of a core in the archetypal filamentary IRDC Gl 1.1 1-0. 12 at few arcsecond resolution to determine its physical and chemical structure. 
The data consist of continuum and line observations covering the C^'^S 2 — » 1 line and the methanol 2^ — » Ij. v, = lines at 3mm 
and the methanol 5^ — > 4j v, = lines at 1mm. Our observations show extended emission in the continuum at 1 and 3 mm. The 
methanol 2j^ — > U v, = emission presents three maxima extending over 1 pc scale (when merged with single-dish short-spacing 
observations); one of the maxima is spatially coincident with the continuum emission. The fitting results show enhanced methanol 
fractional abundance (~3 x 10"^) at the central peak with respect to the other two peaks, where it decreases by about an order of 
magnitude (~4-6 x 10"'). Evidence of extended 4.5 /im emission, "wings" in the CH3OH 2^ — > 1^ spectra, and CH3OH abundance 
enhancement point to the presence of an outflow in the East- West direction. In addition, we find a gradient of ~4 km s"' in the same 
direction, which we interpret as being produced by an outflow(s)-cloud interaction. 
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1. Introduction 



ISM: individual objects (G 11.1 1-0. 12) - ISM: molecules - stars: formation - techniques: interfero- 



The objects now termed Infrared Dark Clouds (IRDCs) were 
first identified by Perault et al. (199^ and |Egan et al. (1998) 
in images of the Galactic plane made with the Infrared Space 
Observatory (ISO) and the Midcourse Space Experiment (MSX) 
satellite, respectively. IRDCs have significant extinctions even 
at 8 and 15 //m and are, thus, seen in silhouette against 
the bright, diffuse, mid-infrared Galactic emission. Millimeter 
and (sub)millimeter molecular studies reveal that some clumps 
within IRDCs have column densities of > 10^^ cm"^, volume 
gas densities of n > 10^ cm~^, and tem peratures T < 20 K 
(e.g.. lCarev et al.lll998l; IPillai et al.ll2006al) . Several studies sug- 
gest that IRDCs have the potential of harboring the earliest 
stages of massive star formation and indeed e vidence for this 
is found toward distinct regions within them (iRathbome et alj 
[2005; Pillai et al. 2006b; Wang et al. 2006). 

IPillai et all (I2006al) mapped the (7, /T) = (1,1) and (2,2) in- 
version transitions of ammonia (NH3) in the IRDC Gl 1 . 1 1-0. 12 
with the Effelsberg 100m telescope and found gas temperatures 
of the order of 15 K for several clumps within this cloud. The left 
panel of Fig. [T] sho ws an overview of the region. Furthermore, 
IPillai et alj (l2006bl) reported the detection of 6.7 GHz class II 
methanol (CH3OH) and 22.2 GHz water (H2O) masers in the 



IRDC Core G 1 1 . 1 1 -0. 1 2P 1 (hereafter G 1 1 . 1 1 P 1 ) . Both CH3 OH 
and H2O masers are known as tracers of massive star formation. 
They found that these mas ers were as sociated with a SCUBA 
dust continuum peak (PI; iCarev et al.l 1200 0) and ascribed the 
kinematics of the methanol masing spots to a maser amplifica- 
tion from a Keplerian disk. 

Recently, Henning et al.l (1201 Ol) observed this filamentary 
IRDC with the PACS (at 70, 100, and 160 /im) and the SPIRE 
(250, 350, and 500 //m) instruments on board of the Herschel 
Space Observatory, with resolutions in the range of ~6"^0". 
For Gl 1.1 IPl, they obtain from a modified blackbody fit to the 
PACS data, a dust temper ature of 24 K, a lu min osity of 1346 
Lq, and a mass of 240 M^^■ lPillai et alj (l2006ah and lLeurini et al.1 
(2007b) have used a two-component model to fit single-dish ob- 
servations of NH3 and CH3OH spectra, respectively, and derived 
temperatures of ~15-18 K for the extended component (with 
size of ~20" corresponding to 0.35 pc) and ~47-60 K for the 
inner component (with a size of ~3" corresponding to 0.05 pc). 

Methanol, which is a slightly asymmetric top molecule, has 
been used to derive p hysical parameters su ch as density and 
temperature in IRD Cs (|Leurini et al.l 2007bl) . high-mass proto- 
stellar objects (e.g.. Leurini et al.T2007al) . massive young stars 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer and the IRAM 30m telescope. IRAM is supported by 
INSU/CNRS (France), MPG (Germany), and IGN (Spain). 
** Member of the International Max Planck Research School (IMPRS) 
for Astronomy and Astrophysics at the Universities of Bonn and 
Cologne. 



(e.g., van der Tak et al. 200 ()) as well as in low-mass protostellar 
system s (iKristensen et alJl2010l) . In particular. Ivan der Tak et al.] 
(I2OOOI) carried out a study toward 13 massive star-forming re- 
gions and found three types of CH3OH abundance (XchjOh = 
NchjOh/Nhi) profiles: XchjOh ~ 10"^ for the coldest sources, 
from lO"** to 10"^ for warmer sources, and 10"^ for hot cores. 
CH3OH has also been associated with outflows where the 
methanol abundance enhancements have been found to be a fac- 
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tor as large as 400-1 OOP (e.g ., lBachiller et al.lll995HGarav et alJ 
l2002HKristensen et al.ll20Tol) . 

In this paper, we present arcsecond resolution millimeter 
continuum and line observations toward Gll.llPl with the aim 
to determining its physical and chemical structure. In Sect. 2, we 
describe our observations carried out with the IRAM Plateau de 
Bure Interferometer In Sect. 3, we present continuum and line 
results together with the analysis. The discussion is presented in 
Sect. 4 and the summary is given in Sect. 5. 



2. Observations and data reduction 

2.1. IRAM PdBI observations 

Gll.llPl was observed with the I RAM six element ar ray 
Plateau de Bure Interferometer (PdBI: lGuilloteau et al.lll992h in 
D and C configurations in 2005 and 2006, respectively. 

During the 2006 observations, one antenna was equipped 
with the prototype New Generation Receiver. Because of a dif- 
ference in the frequency scheme of this receiver, visibilities ob- 
tained in the image sideband in this antenna were discarded. The 
receivers were tuned single side-band at 3 mm and double side- 
band at 1 mm. The 3 mm receivers were centered at 96.64 GHz 
and the 1 mm receivers at 241.81 GHz (see Table [TJ. At 3 mm, 
the C^'^S 2 1 line and the CH3OH 2^ ^ U v, = lines were 
covered by using two correlator units of 80 MHz with a spectral 
resolution of 0.3125 MHz (0.97 kms"') and by two correlator 
units of 320 MHz. The latter were used, free of line, to produce 
the continuum image with a total bandwith of ~450 MHz. At 1 
mm, the CH3OH Sk ^ Ai^ v, - lines were observed with two 
units of 160 MHz with a spectral resolution of 1.250 MHz (1.55 
kms"'). The remaining units of 80, 160, and 320 MHz were 
placed in such a way that a frequency range free of lines could 
be used to measure the continuum flux. The total bandwidth of 
both sidebands was ~700 MHz. 

The phase and amplitude were calibrated with observations 
of the object 1730-130. The bandpass calibration was done with 
3C 273. MWC 349 was used as primary flux calibrator of the 
3 and 1 mm data (see Table [1]). We estimate the final flux den- 
sity accuracy to be ~5% and 10% for the 3mm and 1mm data, 
respectively. Continuum images were subtracted from the line 
data in the visibility plane. The combination of the C and D con- 
figurations provides angular scales from r.'8-12'.'6 (1 mm) and 
4'.'4-32'.'5 (3 mm), i.e., providing information on spatial scales 
of 0.03-0.22 pc and 0.08-0.57 pc, respectively. 

The calibration and data reduction were made within the 
GILDAS packag43 at IRAM Grenoble. Images were created 
with natural weighting and CLEANed using the standard 
Hogbom algorithm. 

2.2. IRAM 30m short-spacing observations 

Gll.llPl was observed in the 2^ v, = (3 mm) and 5a — > 
4k (1 mm) CH3OH bands with the IRAM 30m telescope. An area 
of 90" X 90"was mapped with the SIS BlOO r eceiver at 3 mm 
and 66" x 66" with the HERA receiver (ISchus ter et al. 2004) at 
1 mm with a bandwidth of ~140 MHz. Sampling at 3 mm was of 
15" while at 1 mm was of 6". Conversion from antenna tempera- 
ture to main-beam brightness temperature (Tmb) was performed 
by using a beam efficiency of 0.78 at 3 mm and 0.52 at 1 mm. 
The rms at 1 mm is 0.3 K and at 3 mm is 0.02 K. The data were 



Table 1. Parameters for the IRAM PdBI observations. 



Parameter 



Value 



Pointing center 

Number of antennas 
Baseline range 
Band center 
Primary HPBW 

Synthesised HPBW 

Primary flux density calibrator: 
MWC349 



R.A. (J2000) = 
Decl. (J2000) 



18'U0™28':24 
.i9=22'30'.'5 



24-176 m 
96.64 and 241.81 GHz 
52" at 96.64 GHz 
21" at 241.81 GHz 
6'.'0x3'.'l at 96.64 GHz 
2'.'6x I'.'l at 241.81 GHz 
1.17 Jy at 96.64 GHz 
2.03 Jy at 241.81 GHz 



reduced with the CLASS program, part of the GILDAS software 
package. 

The 3 mm single-dish data are used to create short-spacing 
pseudo-visibilities; these are then merged with the interferomet- 
ric observations. All data are imaged and deconvolved together 
The processing was performed in GILD AS/MAPPING follow- 
ing standard procedures of the UV .SHORT task (see Rodriguez- 
Fernandez et al. 200^ for details on the pseudo-visitbility tech- 
nique). In Sect. l3.3l we present the result of recovering the short- 
spacing information for the 3 mm data. 

Unfortunately, we were not able to merge the interferomet- 
ric and single-dish observations for the 1 mm data due to poor 
signal-to-noise ratio of the latter data. However, we make use of 
the 1 mm single-dish data in the following sections to help in the 
analysis and interpretation of this core. 



3. Results and Analysis 

3. 1 . PdBI: Continuum emission 

3.1 .1 . Morphology and core size 

The right panel of Fig. [1] presents the continuum images at 
1 and 3 mm of Gll.llPl. At the resolution of both images, 
the emission is extended; the peak intensity positions coincide 
with each other and with the maser positions as well. The mor- 
phology at 3 mm resembles that previously imaged with the 
Ber keley-IUino is-Maryland- Association (BIMA) interferometer 
by Pillai et al.1 (i2006bl) with a slightly lower resolution (8'.'3 x 
3'.'9). 

In Table |2] we list the parameters of the continuum emis- 
sion. Deconvolved sizes have been obtained by fitting two- 
dimensional Gaussians to the continuum maps and yielding a 
core size of ~0.16 pc x 0.08 pc with a PA 24° at 3 mm and 
~0.09 pc X 0.07 pc with a PA 165° at 1 mm. We have assumed a 
kinematic distance of 3.6 kpc. 

3.1.2. Mass 

Assuming that the mm continuum is mainly due to optically thin 
dust emission, we estimate the gas m ass, Mg, by using the fol- 
lowing expression (iHildebrandll 19831) 



gd 



(1) 



http://www.iram.fr/IRAMFR/GILDAS 



- http://iram.fr/GENERAL/reports/IRAMjnemo.2008-2-short- 
spacings.pdf 
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18''l0"'40' 20' 18^0'^29' 28' 

Right Ascension (J2000) Right Ascension (J2000) 



Fig. 1. Left: SCUBA 850 yum dCarev et alfcOOOb contours overlaid on the MIPSGAL 24 //m image of the IRDC Gl 1 . 1 1-0. 12. Right: 
PdBI 1.2 mm continuum emission (contours) of Gll.llPl overlaid on the 3.1 mm continuum emission (grey scale image). First 
contour and contour spacing for the 3.1 mm emission are 0.9 mJy beam"' (3cr), the dotted contours show the negative emission 
(-3cr); the synthesized beam (6'.'0 x 3'.' 1 ; PA = 176°) is shown in the bottom left corner First contour and contour spacing for the 1 .2 
mm emission are 6.3 mJy beam"' (3cr), the dashed contours show the negative emission (-3cr); the synthesized beam (2'.'6 x I'.'l; 
PA = 176°) is shown in the bottom right corner The filled triangle at R. A. = 18''10™28'.29, Decl. = -19°22'3a.'5 and the filled squai-e 
at R.A. = 18''10™28!25, Decl. - -19°22'30'.'45 (J2000) indicate the water and (integrated emission) methanol masers, respectively, 
reported by Pillai et al. (2006b). 

Table 2. Parameters of the observed continuum emission. 



Frequency 


Peak Intensity Position 


Peak Intensity 


Flux density_^ 




Deconvolved angular size_^ 


Mass" 


(GHz) 


R.A. (12000) Decl. (12000) 


(mJy beam"') 


(miy) 






{M~ 


96.7 


18''10"'28'.28 -19°22'30'.'5 


7.0 


20.2 


(8'.'9 


t a' 6) X (4'.'4 ± a'5); + 24° ± T 


37 


241.8 


18''10"'28!29 -19°22'3a'5 


30.8 


208.6 


(4'.' 9 + 


a'6) X (3'.'9 ± a'5); + 165° ± 36° 


13 



" Flux density was obtained by summing the continuum emission above 3cr level within a polygon encompassing the core. 
* Major axis x minor axis, at FWHM; position angle of major axis from fits of elliptical Gaussian. 
' Assuming a distance of 3.6 kpc and dust temperature of 60 K. 



where Fy is the observed integrated flux density in Jy, D is the 
distance, R^^ is the gas-to-dust mass ratio, Ky is the dust opacity 
coefficient, and By{Td) is the Planck function at the dust tem- 
perature (Td). The 1 mm flux density was obtained by summing 
the 1 mm continuum emission above 3cr (> 6 mJy beam"') level 
within a polygon encompassing the core; the integrated flux den- 
sity is ~209 mJy. We assume a gas-to-dust mass ratio of 1 00, 
and adopt a /Cy = 1 cm^ g~' (lOssenkopf & Henniiigl ll994t). at 
1mm, for an MRN dMathis et al.lll977l) graphite-silicate grain 
mixture with thick ice mantles, at a gas dens ity of 10^ cm"- ^. 
If we use Td = 60 K ( iLeurini et alj r2007b; Pi Uai et alJ[2006ah . 
a proper temperature for a region where methanol emission is 
arising, a mass of 13 at 1 mm is obtained for Gll.llPl (see 
Table |2|l. This is equiva lent to adopt a dust o pacity coefficient 
ky = 0.1 (v/10'2 Uzf (lBeckwithetalJll990h . which includes 
the gas contribution to the total mass for a R^d = 100, with /3 - 
1.6. Given the uncertainties in the models of dust opacities and 
variation of the gas-to-dust ratio, /3 is within the errors compara- 
ble to the upper limit we obtain in Sect. 13.1.31 

Comparing the flux from the Bolocam Galactic Plane Survey 
(BGPSB Aguirre et al. 2 1 1) 1 . 1 mm with the PdBI 1 .2 nmi con- 



3 http://irsa.ipac.caltech.edu/data/BOLOCAM.GPS/ 



tinuum flux we find that the interferometric observations filter 
out about 75% of the emission. 

To make a rough mass estimate of the core from the 3 mm 
data, all parameters remain the same as with the 1 mm data ex- 
cept for the flux density, ~20 mJy, and Ky = 0.2 cm^ g"' (extrap- 
olating from Table 1 Column 9 in .Ossenkopf & Henning.,19941) . 
We obtain a core mass of 37 Mq. 

3.1.3. Density, temperature, and spectral indeces 

For deeper insight into the source structure, we have analyzed 
the 1 mm and 3 mm continuum data in the u-v domain avoid- 
ing the cleaning and u-v sampling effects. Figure |2] shows the 
averaged amplitudes averaged over a concentric annulus (from 
the continuum peak positions) versus deprojected u-v distance in 
units of wavelength. The amplitudes have been averaged vecto- 
rially and the 3 mm u-v amplitudes in Jy were scaled by a factor 
of 15 to match the 1 mm visibilities. Error bars are Icr statistical 
errors based on the standard deviation of the mean of the data 
points in the bin. 

Since the index of the power law in u-v distance is related 
to the index of the power law in radial distance, r, following 
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10 



10 




1 mm 

3 mm 



Table 3. Molecular parameters of detected line transitions. 



10 



uv Distance (kA) 



10" 



Fig. 2. Continuum emission from Gll.llPl at 3.1 mm (filled 
triangles) and 1.2 mm (filled circles) in the Fourier domain. 
Amplitudes are averaged in bins of deprojected u-v distance 
from the continuum intensity peaks. The 3.1 mm u-v amplitudes 
have been scaled by a factor of 15. The straight solid line rep- 
resents the least-squares fit to the data. Error bars are for 1 cr 
statistical errors based on the standard deviation of the mean of 
the data points in the bin. 



iLoonev et aP (l2003h . we can write F{r) oc v(s) oc 

j(p+9)-3^ where F{r) is the flux density in the image domain, V(s) 
is the Fourier transform of the flux density, s is u-v distance and, 
p and q are the power-law indices of the density (n oc r^P) and 
temperature structure (T oc r"*), respectively. 

The least-squares fit (see Fig.|2]l to both the 1 mm and 3 mm 
data yields a slope of -1.01 + 0.02, corresponding to /:> + <7 ~ 
2.0. Points starting at 40 kA, which corresponds to an angu- 
lar scale of ~6" (0.1 pc), were excluded from the fit because 
the signal-to-noise ratio g ets very low. If the index q = 0.4 
dGoldreich & Kwanll 19741) . the density profile index is 1.6. 

We estimate a mean volume gas density of the core (n - 
where R is the core radius) of 7 x lO*" cm""* using the 



mass derived from the dust observations at 1 mm (13 Mq) and 
the geometric mean of the deconvolved linear size as the core 
radius (R = yl{FWHM^.^^I2) x (FWHM^Jl) = 0.04 pc). 

The scaling factor of 15 corresponds to a spectral index 
a = 3.0 + 0.4 and we thus obtain an opacity spectral index 
yS = 1.0 + 0.4. The derived dust opacity index corresponds, 
taking the upper limit, to a "normal" interstellar dust material 
value. Nonetheless, if taking the lower limit, the dust opacity 
index goes down to a value that is usually found t oward disks 
(e.g., iBeckwifli & Sargent! [19911 iNattaet al.ll2004l) . which has 
been ascribed to grain growth. Alternatively, the presence of 
winds/outflows (see Sect. I4.2l l could mimic a low value of fi 
dBeuther et al.l 120071) . Here we adopt a y8 = 1.6 in the rest of 
the paper. 



3.2. PdBI: Molecular emission 

Spectroscopic parameters of the detected C^'^S 2 1 and 
CH3OH 2k Ik and 5k 4k transitions toward Gll.llPl are 
listed in Table [3] The interferometric integrated intensity maps, 
overlaid on their corresponding continuum emission, along with 
sample spectra are presented in Figs. |3j]5] 



Molecule^ 


Transition 


Frequency 


EJk 








(MHz) 


(K) 


(D^) 




2->l 


96412.9495 


6.94 


7.6678 


CH3OH 


2_i ^ l_i E 


96739.3620 


12.55 


1.2134 




2o ^ lo A 


96741.3750 


6.97 


1.6170 




2o ^ lo £ 


96744.5500 


20.10 


1.6166 






967S'S 'S1 10 


28.03 


1.2442 




5o ^ 4o £ 


241700.2190 


47.96 


4.0402 




5_i ^ 4_i E 


241767.2240 


40.41 


3.8824 




5o -» 4o A 


241791.4310 


34.83 


4.0430 




5±3 ^ 4±3 A 


241832.9100 


84.68 


2.5775 




5i ^4i £■ 


241879.0730 


55.91 


3.9802 




5-2 4_2 E 


241904.1520^ 


60.76 


3.3986 




52 ^ 4, £■ 


241904.6450*" 


57.11 


3.3558 



" Rest frequencies, upper level energies and n~ S from the Cologne 
Datab ase for Molecular Spectroscopy (CDMS; iMuller et all 120011. 
'20051) as of December 2010. 

* Blend of lines. 



3.2.1. C^'^S emission 

The integrated intensity map of the C^"^S 2 — > 1 line, is shown 
in the top panel of Fig. [3] This emission is spatially coincident 
with the 3 mm continuum; moreover, some weaker emission (at 
3cr level) extends along the E-W direction. Observed parame- 
ters toward the peak position in the integrated intensity map are 
presented in Table |4] 

Virial mass. The virial mass of a core with a power-law 
density distribution is Mvir - kicrlR/G (iMacLaren et al.lll988h 
where k[ - (5 - 2p)/{3 - p), cr^ is the three-dimensional root- 
mean-square velocity which is related to the FWHM line width, 
AV, through cr^ - (3/81n2)AV^, R is the core radius and G is the 
Gravitational constant. For p - 1.6 (see Sect. I3.1.3I I. the virial 
mass can be expressed as 



161 X 



AV 

km s" 



Mo 



(2) 



A virial mass of ~135 Mq is calculated when using AV = 
4.1 kms ' (see Table |4]l of the C^'^S optically thin line and R - 
0.05 pc, which corresponds the geometric mean of the decon- 
volved linear size. 

Virial parameter. The virial parameter (iBertoldi & McKed 

for a spherical cloud is defined as = Myii/Mtot, where 



37 M„ 



Mtot is the total mass. If we assume that Mtot 
then ttvir = 3.6. 

It is worth noting that the C^'^S line width could be aff'ected 
by unbound motions, e.g., outflows, partially broadening the line 
and thereby increasing the virial mass estimate significantly and 
the virial parameter as well.^ 

Jeans mass. Following [stabler & PaUal (l2005h . we can cal- 
culate the Jeans mass 



M, 



pl/2 G3/2 



1.0 X — -\ Mo, (3) 

llOK/ \104cm-3j °' ^ ' 



where a is the isothermal sound speed and p is the mass density. 
For n = 7 X 10^ cm"^ (see Sect. lTOl l and T = 60 K, we obtain 
a Jeans mass of ~1.8 Mq. We find a gas mass (Mg) to the Jeans 
mass (My) ratio larger than unity; other i nterferometric studies 
(iRafliborne et al.ll2008l; Izhang et al.ll2009l) toward IRDCs report 
My /My > 1 as well. 
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96.74 96.75 96.76 
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o 

(D 
Q 



o 

CM 



o 

CM 
CM 



O 

CM 
CM 



O 

CM 
C^ 



CH3OH 2i-liE 








18 10 30 28 

Riqht Ascension (J2000) 




96.74 96.75 96. 7f 
Rest Frequency (GHz) 



(3.7,-7.3) 




96.74 96.75 96.76 
Rest Frequency (GHz) 



Fig. 4. In grey-scale the 3.1 mm continuum emission of Gil . 1 IPl . Top middle: Contour image of the CH3OH emission integrated 
under the line 2i liE. First contour and contour spacing are 0.66 Jy beam"' kms"' (3cr). The synthesized beam (6'.'0 x 3'.'2; 
PA - 174°) is shown in the lower right comer Bottom middle: Contour image of the CH3OH emission integrated under the 2_i — > 
l-\E, 2o — > IqA, and 2o — > IqE transition lines. First contour and contour spacing are 0.45 Jy beam"' kms"' (3cr). The dashed 
contours show the negative emission (-3cr). Shown are also spectra at four different positions; positions are given in each inset in 
parenthesis in units of arcsec. 



Table 4. Observed molecular line parameters from Gaussian fits. 



Molecule 


Transition 


vlsr 


Peak Tub 


AV (FWHM) 








(kms-') 


(K) 


(kms"') 


(Kkms"') 






Offset (0",0") 








2^ 1 


29.90 (0.04) 


1.31 (0.04) 


4.1 (0.1) 


5.6 (0.1) 






Offset (l'.'5,0") 






CH3OH 


2_i l_i E 


29.85 (0.05) 


1.49 (0.04) 


5.1 (0.2) 


8.0 (0.2) 




2o ^ lo A 


29.72 (0.04) 


1.89 (0.04) 


3.6 (0.1) 


7.3 (0.2) 




2o ^ lo £ 


29.87 (0.09) 


0.78 (0.04) 


4.4 (0.2) 


3.6 (0.1) 




2, li £ 


29.79 (0.01) 


0.49 (0.04) 


3.1 (0.2) 


1.6 (0.1) 




5o ^ 4o £ 


29.71 (0.20) 


0.70 (0.04) 


4.1 (0.5) 


3.1 (0.3) 




5_i ^ 4_i £ 


29.92 (0.09) 


1.48 (0.04) 


4.1 (0.2) 


6.5 (0.3) 




5o ^ 4o A 


30.00 (0.05) 


1.80 (0.04) 


4.3 (0.2) 


8.3 (0.3) 




5i ^ 4i £ 


29.78 (0.41) 


0.33 (0.04) 


3.9(1.3) 


1.4 (0.3) 



Notes. CH3OH 2<, —> III and 5^. — > 4^. line parameters listed in this table were used in the rotational diagram method described in the text. The 
CH3OH 5k — » 4k transitions have been convolved to the 2j. — ^ U resolution. 



From our interferometric observations at 1mm, and taking a 
flux density equal to the 3cr detection level, we derive a mass 
sensitivity limit of 0.4 Mq. Comparing this value with the Jeans 
mass, we note that our mass sensitivity limit is good enough to 
detect fragments of the order of the Jeans mass. 



3.2.2. CH3OH emission 

Figure |4] shows the integrated intensity map of the CH3OH 
2jt — > U quartet of lines (bottom middle). This emission presents 
three maxima; one of them is spatially coincident with the 1 
and 3 mm continuum emission which peaks at (l'.'5,0"), the 
other two are located at ( 1 3'.' 1 , 0") and (-9'.'5, 2'.'2). We plot sam- 
ple spectra at four different positions, including the three max- 
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Fig. 3. Top: Contour image of the emission integrated un- 
der the C^"^S 2—> 1 line, in the velocity range from 22.4 to 
37.9 kms"' , overlaid on the grey-scale 3.1 mm continuum 
emission of Gll.llPl. First contour and contour spacing are 
0.09 Jy beam"' kms"' (3cr), the dashed contours show the neg- 
ative emission (-3cr). The synthesized beam (6'.'0 x 3'.'2; PA - 
176°) is shown in the lower right corner Bottom: Spectrum of 
C^'^S 2— > 1 toward the peak position on the integrated intensity 
map: R.A. = 18'n0'"28!24, Decl. = -19°22'30'.'5 (J2000). 



ima and one position on the envelope (3'.'7, -7'.'3), in order to 
compare line profiles, absolute line strengths, and the relative 
strengths among lines. The high excitation CH3OH 2i — > liE 
transition that arises from a smaller region is also shown (top 
middle). 

Line widths at and close to the three peak positions are, 
in average ~5 kms"' with line profiles showing red- and/or 
blueshifted "wings", while in the outer parts of the envelope, 
e.g., at the (3'.'7, -7'.'3) position, the line width is narrower (~2 
kms"'; see spectra inFig.|4|i. 

In Fig. |5] we present the integrated intensity map of the 
strongest line CH3OH 5o — > A-qA. The emission distribution is as 
extended as that of the 1 mm continuum and the CH3OH 2i —> 
liE emissions. 

Local thermodynamic equilibrium (LTE) analysis. We 

derived the rotational temperatures and methano l column densi- 
ties by using the rotational diagram method (e.g. .'Cummi ns et aP 
[l986; Goldsmith & Langer 1999). Briefly, for an optically fliin 
line, the column density in the upper level (A^„) can be written as 
Nu - 3kW/{^n^vS p^), where k is the Boltzmann constant, W is 
the integrated intensity of a line ( J TMsdV), v the frequency, S 
the line strength, p the dipole moment. Under LTE conditions, 
the population of each level follows a Boltzmann law at the gas 
temperature T hy Nu - (N/Q)giie^^"^''^ , where is the total 
column density, Q the partition function, and g,, and £„ are the 
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Fig. 5. 7b;?: Contour image of the emission integrated under the 
5o 4oA transition line, overlaid on the grey-scale 1 .2 mm con- 
tinuum emission of Gll.llPl. First contour and contour spac- 
ing are 0.16 Jy beam"' kms"' (3cr); the dashed contours show 
the negative emission (-3cr). The synthesized beam (2'.'6 x I'.'l; 
PA - 174°) is shown in the lower right comer. The dotted cir- 
cle indicates the interferometer primary beam (21") at this fre- 
quency. Bottom: Spectrum of CH3OH Sk —> 4-k towards the po- 
sition: R.A. = 18''10™28H0, Decl. = -19°22'31'.'9 (J2000). The 
HNCO llo.ii -> lOo.io line (at 241774.0320 MHz with EJk = 
69.67 K and p^S = 27.4590 D") is shown as well. The "G" shows 
an artifact produced in the central channels of a correlator unit 
due to Gibbs effect. 

statistical weight and the energy of the upper level, respectively. 
We then obtain 



log 



log 



N\ loge E„ 

Q 



T k 



(4) 



Using Eq. |4l we plot log(Nu/gii) versus E,i/k (see Fig. |6]l 
and make a least-squares fit for T and taking into account 
both CH3OH 2i —> It and 5^ —> 4^ bands. The partition func- 
tion of CH3OH is e = 2[pi^]'''2 ( |Turnedll991h . where h is 
the Plan ck constant and A, B, and C are the rotation constants 
(CDMS: iMuiler et al.ll200ll l2005h . A g = 1.23 r'-^ was used 
in the calculations. The CH3OH 5k — » 4^ data were smoothed to 
the resolution of the —> Ik data. CH3OH transition lines that 
were used in this fit are listed in Table |4] 

To estimate the CH3OH fractional abundance (relative to 
H2), ^cHsOH, we calculate the beam-averaged H2 column den- 
sity, A'h, , following the expression 



.peak n 



Ky By(Td) O.PH2 mn 

^peak 



(5) 



where /j; is the pe ak intensity in Jy/beam, p^^ - 2.8 
dKauff^mann etalJl2008h is the molecular weight per hydrogen 



6 



L. Gomez et al.: Observations of methanol in the IRDC core Gll.l 1-0.12P1 





m 




CH3OH 2k- 




13 








2.5 


- 






- 


12 










1 .5 


- T,„,= 15K 


15 -2 

cm 




1- 



■10 



Fig. 6. Rotational diagram of CH3OH (2-1) and (5^) transition 
lines towards the (r.'5,0") position. The straight line is a least- 
squares fit to the data corresponding to the rotational temperature 
and column densities indicated in the bottom left. Error bars are 
for Icr uncertainties. 

molecule, and Q = (7r6',ni„0maj)/(41n2) is the beam solid angle of 
an elliptical Gaussian beam with minor and major axes flmi,, and 
^maj, respectively. 

We obtained a rotational temperature of 15 K and a CH3OH 
column density, A^chjOh, of 1.2 x 10'^ cm"^ towai'ds the (1'.'5, 0") 
peak position. Since the lines are likely sub-thermally excited, 
we take a dust temperature = 60 K (see LVG analysis) and 
obtain a A^h^ = 4.6 x 10^^ cm"^ (from our 1 mm continuum data 
convolved to the 3 mm data) and a Xch,oh of 2.6 x 10"^. 

At the (13'.'1,0") position, we get a Xch,oh = 3.7 x lO"'', 
when using A^ch,oh = 2.9 x 10'"* cm^^ and A^h, = 7.7 x lO^^ 
cm"^. A rotational temperature of 13 K was calculated. At the 
(-9'.'5, 2'.'2) position, the values for the Xch,oh, A^chjOh, and 
A^H, are: 5.6 x 10"^ 3.8 x 10'^ cm"2, and 6.8 x 10^^ cm^^, 
respectively. For the (-9'.'5, 2'.'2) peak, we fixed T = 15 K be- 
cause of low signal-to-noise of one of the three lines used in the 
fit. A^CHiOH for both peaks (13'.' 1,0") and (-9'.'5, 2'.'2) were cal- 
culated from the 3 mm interferometric map only and A^h, from 
the SCUBA 850 fj.m dust continuum; the interferometric obser- 
vations were convolved to the SCUBA beam size ( 14"). 

Large velocity gradient (LVG) analysis. iLeurini et"al] 
(12004 ) found several ratios, from some CHsOH-E transitions, to 
be calibration-independent tracers of density in the 1 and 3 mm 
bands. We have used their LVG calculations and have plotted the 
following integrated intensity line ratios: 5() — > 4o/5_i — > 4_i, 
5i 4i/5_i — > 4_i, and 5i — > 4i/5o — > 4o, in a temperature- 
density (T-n) plane. In Fig. |7l the observed ratios, towards the 
(l'.'5,0") position, are represented by the black dash-dotted line, 
blue solid line, and pink dashed line. The 2^ 1^ quartet of 
lines in the 3 mm band were not used in the analysis because 
they are blended. 

Given the results for CH3OH column densities in the 
LTE analysis, the explored LVG models for the excitation of 
methanol run for temperatures in the range of 10-200 K, for 
densities 10^-10*^ cm"^^, and for two CH3OH column densities 
per line width lO'** and lO'^ cm ^/(kms '). 

For an average line width of 4 kms \ this corresponds to 
A^CHjOH-E = 4 X 10'"* cm"^ and 4 xlO^^ cm"^, respectively. 
Depending on the column density, two possible solutions are 
found: at A^chjOh-e = 4 x 10'** cm"^, the three integrated in- 
tensity line ratios intercept at T = 17 K and « = 10^ cm"^. 
On the other hand, at A^ch^oh-e = 4 x 10'^ cm"^, the inter- 
ception is at r = 60 K and « = 2 x 10^ cm""* (see Fig. |7]|. 
The methanol column density inferred in this case is consistent 
with the LTE analysis within a factor of 3, but the intensities in 
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Fig. 7. Results of statistical equilibrium calculations for 
CH30H-£'. The 5o ^ 4()/5_i 4_i (black dash-dotted 
line), 5i —> 4i/5_i —> 4_i (blue solid line), and 5i —> 
4i/5o — » 4o (pink dashed line) observed integrated inten- 
sity line ratios towards the (l'.'5,0") position in logarithmic 
scale are shown as function of H2 density and temperature at 
A^ch,oh/AV = 1014 cm-2/(kms-i) (top) and A?ch,oh/AV = lO'^ 
cm"^/(kms"') (bottom). 



this model are one order of magnitude higher than the observed 
line intensities. In other words, to fit the observations with the 
A^ch,oh-e = 4 X lO'*" cm"^ model, we need a beam filling factor 
of 0.1 resulting in a more compact source (with a size of 0'.'5 
corresponding to 0.009 pc or 1800 AU) , even smaller than th e 
one modelled bv Pillai et al.l ( l2006al) and lLeurini et al. I (l2007bh . 

To gain more insight into the possible T-n degeneracy, in 
Figure |8] we show n vs. A^chjOh plots, at different temperatures 
(10-70 K), where the color-coding for the ratios (along with their 
Icr values) are the same as in Fig. |7j additionally, we show in 
each panel the peak main beam brightness temperature of the 
5o — » '^oE transition (green dotted line). The parameters ob- 
tained in the regions in the left side of the green dotted line imply 
beam filling factors greater than 1, which is not possible. In each 
panel, the areas of interest are those that lie in the right side of 
the green line and, at the same time, are delimited by the inte- 
grated intensity ratios (taking into account their Icr values). We 
see that, for T in the range 30-70 K, the n remains constant up 
to NcH.OH ~ 4 X 10"^ cm.-^, while for the very low T (= 10 K), 
there is no solution. 

Given our density estimate from the dust continuum, we 
view the low T, high n solutions with a large filling factor as 
unlikely. Therefore, although we cannot strongly constrain the 
physical parameters of the gas emitting methanol from the LVG 
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Fig. 8. Results of statistical equilibrium calculations for CHjOH-is. Observed integrated intensity line ratios are color-coded as 
in Fig. Q and plotted in logarithmic scale (along with their Icr values) as function of H2 density and A^ch,oh/AV at different 
temperatures (10-70 K). The green dotted line represents the peak main brightness temperature of the 5o — » 4-()E transition. 



analysis because of the limited number of observed transitions, 
we can exclude the low temperature solution and put a lower 
limit to the kinetic temperature of the gas to 30 K. This result 
is in agreement with the detection of methanol maser from this 
region and with a previous measurement of T = 60 K from NH3 
observations. Therefore, we adopt T = 60 K as kinetic tempera- 
ture at the peak of the continuum emission. 

3.3. PdBI+30m: extended CH3OH emission 

A comparison of the interferometric data with the single-dish ob- 
servations shows that only ~30% of the CH3OH 2^ — > 1^ mea- 
sured with the 30 m telescope is imaged in the PdBI observa- 
tions. The combined PdBI and 30 m telescope data recover al- 
most 100 % of the flux at the central (l'.'5,0") position. In Fig.|9l 
we show the resultant map after combining the CH3OH 2^ — » 
lie single-dish observations with the corresponding PdBI data. 
The methanol emission is extended already in the interferomet- 
ric observations (mapping angular scales as large as 32'.'5) and 
confirmed to be so in the resultant merged map with the single- 
dish data. It is important to note that the emission is extending 
over 1 pc scale in the merged map (see Fig.|9). Similarly, silicon 
monoxide (SiO) was found to be widespread oyer a 2 pc scale in 
the IRDC G035.39-00.33 (Ji menez-Seiraet al.ll2010l) . The 2i ^ 
liE transition was not detected in the single-dish observations, 
which is expected since this emission does not appear to be ex- 
tended in our interferometric data. 

As for the 1mm data, the missing flux in the interferometric 
observations is ~70%. In the 1 mm single-dish observations, two 
of the strongest transitions showed up, namely: the 5_i — > 4_i£' 
and the 5o — > 4oA lines. Both lines are above 3cr noise level at 6 
positions out of 144 observed positions. 

How is the derivation of column densities and rotational tem- 
peratures affected by missing flux? Using the combined data at 



3 mm only, the slope, i.e., the rotational temperature, in the rota- 
tional diagram is not affected, whithin the errors, while CH3OH 
column densities are higher by a factor of two. 

4. Discussion 

4. 1 . Properties ofG11.11P1 compared to other cores 

Despite the fact that our interferometric observations have 
enough sensitivity to resolve the Jeans mass (1.8 Mq) at 1 mm, 
the angular resolution (2'.'6xl'.'l) is smaller than but very close 
to the Jeans length of 3'.'2 (0.05 pc). Therefore, it is difficult to 
conclude whether Gl 1.1 IPl does fragment or not into conden- 
sations (referred as substructures within the cores) as is often 
found when zoo ming into these cores (e.g.. lZhanget all 120091: 
iBontemps et alJlIOlQ) . If Gll.llPl does not fragment would 
mean that it has already accumulated most of its mass that will 
go into the final star even though ~75 % of the (dust) emission is 
filtered out by the interferometer With higher angular resolution 
observations, it will be possible to support either pos sibilities. 

The virial parameter of Gll.llPl obtained bv iPillai et al.l 
(I2006ah . from single-dish NH3 and dust continuum (at 870 //m) 
observations, is unity on the scale of ~0.9 pc. Although, in the 
present work, we find that Mvh is much larger than the Mg on 
the scale of ~0. 1 pc which would indicat e that gravity is not the 
dominant force (lBertoldi(S:McKeell992h . We caution, however, 
that the high value of the virial parameter (avir - 3.6) could be 
due to the missing flux in the interferometric data and to the 
fact that the virial mass is estimated from the C^'^S line width, 
whose broadening i s likely in part produced by outflows. In fact, 
iPillai et aP (l2006lj) reported an increase in line widths, from 
1.29 kms"' in the NH3 (1,1) line to 3.2 kms"' in the NH3 (3,3) 
data. The increase in NH3 line width toward the higher (J, K) 
transitions suggest a strong influence of star forming activity on 
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Fig. 9. CH3OH 2ii —> Ik V, = Q integrated intensity map com- 
bining PdBI and 30m data. The emission is integrated over the 
2_i l-iE, 2o — » IqA, and 2o — > IqE transitions. First con- 
tour and contour spacing are 0.39 Jy beam"' km s"' (3cr), the 
dashed contours show the negative emission (-3cr). The synthe- 
sized beam (6'.'2 x 3'.'3; PA - 176°) is shown in the upper right 
corner. The dotted circle indicates the interferometer primary 
beam (52"). 



the line width at smaller scales. The broad line width in the C^^^S 
from PdBI is consistent with that notion. 

Based on the SCUBA 870 //m column density, we estimated 
a surface density of S = 0.4 g cm"^ for G 11 . 11 P 1 , which is close 
to the thresholds of 0.7-1.5 g cm"^ (or 3 320-7 1 10 Mp pc'^) re - 
quired to forni stars of 10-200 Mq (iKrumholz & McKedllOOg) . 
It is important to point out that these authors were interested in 
clouds where no massive stars have yet formed while Gl 1 . 1 IPl 
shows already signposts of star formation activity. 

We have calculated a density profile power-law index of 
1.6 from interferometric observations at 1 and 3 mm; sim- 
ilar values, in the range between 1.5 and 2, for cores in 
the IRDC G28.34-h0.06, have been reported by Zhang et al.' 
(I2009.) and for embe dded protostellar so urces in the protocluster 
IRAS 05358H-3543 (iBeutheret al.ll2007l) as well. 



4.2. The disk/outflow system? 



iPillai et alj (l2006bl) found a velocity gradient of the CH3OH 
maser emission and explained it as the maser spots being located 
in a Keplerian disk. Interestingly, the spread of the maser com- 
ponents is in the North-South direction which is perpendicular to 
the CH3OH 2ii — > Ik emission that we propose is originated by 
an outflow (in the East-West direction). The spread of the masing 
spots is translated into a radius of ~450 AU for a hypotetical disk 
inGll.llPl. Note that candidate disks in high-mass protostars 
have masses in the range 0.2^0 Mq and radii of 500-10000 
AU (ICesaroni et al.ll2006l) . Our interferometric data set is capa- 
ble of imaging the biggest disk structures, but not those inferred 
in|Pillai et al. (2006b). 

Gll.llPl has been cataloged bv ICvganowski et"an (120081) 
as a "possible" massive young stellar object (YSO) outflow can- 
didate. This categorization was based on the angular extent of 
the extended excess 4.5 fim emission, i.e., the extent of green 
emission in a three-color RGB image. 

More evidence pointing to the presence of outflows comes 
from the non-Gaussian methanol line profiles (see Fig. |4] and 
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Fig. 10. Position-velocity plot of the methanol 2k Ik lines, 
cutting along the PA=96° on the map shown in the middle bottom 
panel of Fig.H] Shown are the transitions: 2_i — > \-\E, 2o — > 
IqA, and 2i 1 \E. The offset is measured positive towards East 
and from the pointing center of the observations. First contour is 
0.018 in steps of 0.036 Jy beam"', the dashed contours show the 
negative emission (-0.018 Jy beam"'). 



i Leurini et al.l2007"bl) . In Fig.[TO] we analyze the molecular emis- 
sion in a position-velocity plot and note that there are "wings" 
of redshifted emission towards negative off'sets and blueshifted 
emission towards positive offsets (see also Fig. |4|l. Moreover, 
we find a gradient of ~A kms"', in the east-west direction, for 
the CH3OH 2k — » 1a emission and speculate, it is due to an 
outflow(s)-cloud interaction. The C'''*S 2— >1 and CH3OH 5^ — > 
Ak lines are detected mainly towards the central peak/dust emis- 
sion, and no velocity gradient is found. 

Another possibility is that the CH3OH 2k ~> U emission 
is actually tracing a toroid (see ICesaroni et al ] l2007h like in the 
interferometric observations of CH3O H 60 — > 5o towards the 
IRDC 18223-3 dFaflscheer et alj2009l) . However, the gradient of 
the maser spots is orthogonal and not aligned with the gradient 
of the methanol emission we present in Fig. [TO] 

4.3. CH^OH abundances 

Due t o the fact that methano l molecules can be sub-thermally ex- 
cited dBachiller et al.|[T995h . flie temperatures derived from the 
rotational diagram are low (~15 K) even at the central peak 
where masers have been detected. To overcome this issue we 
have used the LVG method and found that a higher temperature 
is more reliable. 

Purely gas-phase model calculations p redict Xch,oh to be of 
the order 10"'^-10"'° (iGarrod et al.' 2006), which cannot repro- 
duce the high values we report in this paper Gan^od et al. (20()^ 
conclude that the production of methanol is carried out on the 
surfaces of dust grains followed by desorption into the gas. Then 
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the XCH3OH values in Gil. P 11 could be produced by desorption 
of icy mantles on the dust grains. 

As found in Sect. 13.2.21 the highest Xch,oh in GU.UPl is 
at the central position (XchjOh ~ 3 x 10"^), and it decreases 
by more than one order of magnitude at the other two peaks, 
where the Xch,oh (~ 4 - 6 x 10"^); the s ame trend was foun d 
also by LVG-modeling of single-dish data dLeurini et al .I2007bl) . 
Our abundances are enha nced compared to d ark clouds values 
(^CHjOH ~ 10 - 10'^: iFriberg et al.lll988h . also supporting 
the idea that methanol in Gll.llPl is forming mainly through 
non-thermal desorption by the presence of outflow(s). 

5. Summary 

We have performed continuum and line observations with the 
PdBl towards the IRDC core Gll.llPl. Our main results can be 
summarized as follows: 

- The analysis of the mm continuum maps provides a very de- 
tailed physical structure of the core. 

- Evidence of extended 4.5 fim emission, "wings" in the 
CH3OH 2i — > lit spectra, and CH3OH abundance enhance- 
ment point to the presence of an outflow in the East- West 
direction. 

- We find a gradient of ~4 kms ' for the CH3OH 2^ U 
emission, which we interpret as being produced by an 
outflow(s)-cloud interaction. 

- The fitting results show enhanced methanol fractional abun- 
dance (~3 X 10"**) at the central peak with respect to the 
other two maxima where the methanol abundance is lower 
by about an order of magnitude (~4-6 x 10"^). 
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